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Aclaración inicial 


El polemista expone errores y amplía conceptos 
teóricos utilizados por el autor, con el objetivo 
básico de dar al lector y al usuario potencial del 
artículo una panorámica completa de los análi- 
sis regionales de frecuencia de gastos máximos 
anuales y del alcance real del estudio. Al autor, 
las observaciones del polemista le pueden ser 
útiles cuando actualice el estudio cada cinco años, 
como lo propone al final de su artículo. "Tales 
observaciones se presentan según tres tópicos: 
(a) aspectos operativos del análisis regional; (b) 
deducción de ecuaciones empíricas del gasto 
medio anual, y (c) estimación del tiempo de 
concentración de las cuencas. 


Errores de omisión respecto a los análisis 
regionales 


Relativo a la calidad estadística de los datos 
hidrométricos 


El polemista destaca que no se expuso. ¿Cómo 
se verificó la homogeneidad, independencia 
y ausencia de componentes determinísticas 
en los registros hidrométricos originales? En 
concreto, ¿cómo se verificó que los 309 regis- 
tros de datos son estacionarios? Lo anterior 
es fundamental antes de proceder a verificar 
la homogeneidad hidrológica regional, ya que 
si existe tendencia en los registros, debe ser 
tomada en cuenta en el análisis regional, como 
lo han mostrado Cunderlik y Burn (2003), y 
Leclerc y Ouarda (2007). 


El polemista indica que como mínimo se 
debió aplicar la prueba de discordancias (Hos- 
king 6 Wallis, 1997; Campos-Aranda, 2010) a los 
registros de datos hidrométricos de cada región 
O grupo de regiones adoptado, para detectar 
registros anómalos y eliminarlos; o bien, para 
comprobar que eran adecuados y continuar con 
la verificación de la homogeneidad hidrológica 
regional. 


Relativo al agrupamiento de regiones y 
estaciones hidrométricas 


En el primer párrafo de la página 22, el autor 
menciona: “Los grupos de regiones hidrológicas 
se establecieron como se indica en el cuadro 1 
(ver figura 2).” 

El polemista formula la pregunta siguiente: 
¿tal agrupamiento se debe a escasez de registros 
hidrométricos o a similitud de comportamiento 
hidrológico? Cualquiera que sea la respuesta, se 
debió indicar cómo fue verificado que tal agru- 
pamiento resultaba conveniente y /o aceptable. 

El polemista indica que el inciso de la página 
18 titulado “Agrupación de regiones”, se debió 
llamar “Agrupación de estaciones hidrométri- 
cas”, ya que ahí se describe cómo se formaron 
los grupos G1, G2 y G3, y cómo se verificaron 
espacialmente en cada región. Esta subdivisión, 
el polemista la considera inadecuada o no con- 
veniente, ya que al formar 18 regiones (cuadro 1) 
y haber procesado 309 estaciones hidrométricas, 
se obtienen en promedio 17 estaciones en cada 
zona geográfica y tal cantidad, afortunadamen- 
te, está próxima al número máximo que se ha 
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recomendado de 20 estaciones por región 
homogénea (Hosking éz Wallis, 1997). 


Relativo a la homogeneidad hidrológica 
regional con datos hidrométricos 


Con excepción del texto del último párrafo de 
la “Introducción”, en el cual el autor indica 
que: “Además, tomando en cuenta la similitud 
en las características estadísticas de los gastos 
máximos registrados en grupos de cuencas, se 
desarrolló un análisis regional que permite....”, 
el polemista destaca que no se verificó la homo- 
geneidad hidrológica de las regiones ni de los 
grupos que con ellas se formaron. Por lo tanto, 
todos los resultados del cuadro 4 del autor son 
cuestionables. Para realizar tal análisis con base 
en los datos disponibles de gasto máximo anual 
se recomienda consultar a Hosking y Wallis 
(1997), Escalante-Sandoval y Reyes-Chávez 
(2002), y Campos-Aranda (2010), además de 
todas las referencias citadas por el autor en el 
cuerpo central de la “Introducción”. 


Relativo a la homogeneidad hidrológica 
regional con fechas de crecientes 


El autor indica en la “Introducción” que el agru- 
pamiento de cuencas se ha basado principal- 
mente en dos tipos de variables para establecer 
la similitud entre cuencas: (1) sus características 
fisiográficas, y (2) las propiedades estadísticas 
del registro de crecientes y de los datos clima- 
tológicos disponibles. 

El polemista ha encontrado (Burn, 1997; 
Cunderlik € Burn, 2002) que el principal ar- 
gumento contrario al uso de las características 
fisiográficas es que no hay garantía de que la 
semejanza fisiográfica implique similitud en la 
respuesta hidrológica, sobre todo si las cuencas 
están distantes, o tienen coberturas de suelo y 
vegetación diferentes. Por otra parte, el pro- 
blema básico con respecto al uso de los datos 
hidrométricos implica una situación espuria 
(Clarke, 2010), pues tales datos se usan para 
formar las regiones y también para verificar su 
homogeneidad. 
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Por lo anterior, el polemista sugiere tomar 
en cuenta una verificación de la homogeneidad 
hidrológica regional con base en la fecha de ocu- 
rrencia media y la regularidad estacional de los 
gastos máximos anuales, ya que una similitud 
en tales parámetros implica una semejanza en 
características fisiográficas y meteorológicas. En 
México, el estudio de la fechas de ocurrencia 
de las crecientes fue introducido por Ramírez, 
Gutiérrez-López y Ruiz-Silva (2009), y ha sido 
aplicado por Campos-Aranda (2014). 


Relativo al análisis probabilístico de las 
series normalizadas 


En opinión del polemista faltó indicar: 


1. Respecto a las funciones de distribución de 
probabilidades (FDP), se debió citar: ¿cuáles 
fueron las utilizadas? y ¿se aplicaron las 
FDP establecidas bajo precepto, es decir: 
Log-Pearson tipo III, General de Valores 
Extremos y Logística Generalizada? 

2. ¿La estimación de los parámetros de ajuste 
de las FDP se realizó con el método de 
momentos, con el de máxima verosimilitud, 
con el de momentos L o con cuál? 

3. ¿Cómo se seleccionó el mejor ajuste de 
las FDP? ¿Se aplicaron pruebas o test 
estadísticos o se empleó el error estándar 
de ajuste? 


Para dar respuesta a las preguntas de los 
incisos 1, 2 y 3, el polemista sugiere consultar 
las referencias: Hosking y Wallis (1997); Rao y 
Hamed (2000); Asquith (2011), y Meylan, Fabre 
y Musy (2012). 


Errores de concepto al estimar MOMIA 


En opinión del polemista, existe un error 
conceptual en seleccionar al tiempo de concen- 
tración (Tc) de cada cuenca como la segunda 
variable predictiva en la estimación de la cre- 
ciente O gasto medio anual (Qma), designado 
por el autor MOMIA. Este error se observa 
en los resultados mostrados por el autor en el 
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cuadro 5, cuyas ecuaciones empíricas en función 
del área de cuenca (4) y del Tc en las regiones 
hidrológicas 25 y 28-29 tienen exponente cero 
en el Tc, y en las regiones 20-21, 26VM, 27 y 30, 
el exponente del Tc es positivo. 

El polemista también considera otro error 
conceptual la selección como cuarta variable 
predictiva de la retención potencial máxima 
(S), pues tal cantidad se estima (ecuación (2)) 
con base en el número N de la curva de escurri- 
miento y tal parámetro hidrológico conlleva una 
estimación subjetiva (Mockus, 1972; Campos- 
Aranda, 2015). 

El autor señala, en el primer párrafo de las 
“Conclusiones y recomendaciones”, que el Qma 
es el concepto eje del método expuesto, ”...de 
tal manera que al relacionarla con características 
medibles en cualquier cuenca (área, tiempo de 
concentración, promedio espacial de la media 
de precipitaciones máximas anuales, retención 
potencial máxima) es posible estimarla de forma 
sencilla, aun para cuencas no aforadas y,...”. 

El polemista opina que únicamente el área 
de cuenca puede ser medida y el resto de va- 
riables predictivas utilizadas por el autor debe 
ser estimada; para el caso del Tc, con base en 
varias ecuaciones empíricas aplicadas sin ser 
extrapoladas, como se explica posteriormente. 

El autor indica en el cuarto párrafo de las 
“Conclusiones y recomendaciones” que *“ 
varios grupos de regiones tuvieron exponentes 
positivos asociados con la retención potencial 
máxima y/o tiempo de concentración, lo cual 
no parece lógico, por lo que se recomienda...”. 

En opinión del polemista, la obtención de 
exponentes positivos en las dos variables pre- 
dictivas citadas (S y Tc) implica una relación no 
funcional o incorrecta de la ecuación empírica 
del Qma, ya que tanto S y Tc, al ser mayores, 
originan un Qma menor y viceversa; por lo 
tanto, deben aparecer como denominadores en 
la fórmula empírica deducida. El autor obtiene 
(cuadro 5) seis ecuaciones empíricas incorrectas 
en el Tc y nueve en S, de las veinte que presenta 
de cada una. 

En resumen, el polemista propone que 
en lugar de utilizar una estimación del Tc se 


debieron emplear características fisiográficas 
de la cuenca, que son estimables por medición 
y que se sabe que de ellas depende el Tc (ver 
ecuación (1) siguiente), como son: la longitud 
del cauce principal, su pendiente promedio, la 
densidad de drenaje, etcétera. Una búsqueda 
de variables predictivas y su selección para 
evitar la multicolinealidad se puede consultar 
en Campos-Aranda (2013). 

El polemista considera otro error de concep- 
to la selección de la tercera variable predictiva, 
denominada volumen de precipitación en la 
cuenca (V), ya que al ser el producto de la media 
de la precipitación máxima diaria (PMD) anual 
por el área de cuenca, la ecuación (5) del autor 
incurre en una correlación espuria (Benson, 
1965). Por otra parte, el polemista indica que 
siendo el Qma una predicción de periodo de 
retorno dos años, es acertado usar, como lo hace 
el autor, la media de la PMD anual como va- 
riable predictiva, pero sólo ella, lo cual ha sido 
verificado en otros países (Prosdocimi, Kjeldsen 
éz Svensson, 2014). 


Error de omisión de los índices de 
desempeño de las ecuaciones empíricas 


En la figura 3, el autor expone para la Región 
Hidrológica núm. 18 una ecuación empírica 
que estima el MOMIA en función del área 
de cuenca, obtenida con 15 datos y con un 
R?= 0.8173. Para esta misma región, en la figura 
5 expone otra ecuación similar, pero obtenida 
con nueve datos y con un R? = 0.7834. El po- 
lemista pregunta: ¿por qué se eliminaron seis 
cuencas de tal región? Este tipo de decisiones 
debe ser explicado y como mínimo indicar el 
intervalo de aplicación de cada ecuación empíri- 
ca deducida (cuadro 5), así como sus índices de 
desempeño, evaluados en el dominio real (ver 
Campos-Aranda, 2013). 


Error operativo al estimar el Tc 
El polemista considera un error operativo esti- 


mar el tiempo de concentración (Tc) utilizando 
una sola ecuación empírica, en este caso la de 
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Kirpich, quizás la menos indicada de las aplica- 
das comúnmente. Lo anterior se intenta aclarar 
en los incisos siguientes, no teniendo como obje- 
tivo presentar fórmulas sino únicamente indicar 
dónde se pueden consultar y qué intervalo de 
aplicación tienen. 


La incertidumbre inherente del Te 


Después del área de la cuenca, el Tc es el pará- 
metro principal de la mayoría de los métodos 
hidrológicos de estimación de gastos máximos 
o bien se utiliza para cuantificar el tiempo al 
pico (Tp) del hidrograma de respuesta o el tiem- 
po de retraso (1,) de la cuenca. A pesar de tal 
universalidad, el Tc, como ha señalado McCuen 
(2009), no tiene una definición computacional 
única, y por ello recabó y expuso seis, y sugirió 
la siguiente como la menos incierta para tal 
estimación: “diferencia entre los tiempos en 
que ocurren los centros de masa de la lluvia en 
exceso y del escurrimiento directo”. Grimaldi, 
Petroselli, Tauro y Porfiri (2012) citan dos de- 
finiciones teóricas del Tc, adicionales a las seis 
operativas de McCuen (2009), las cuales son: 
(a) tiempo que le toma a una gota de lluvia en 
llegar a la salida de la cuenca, partiendo de su 
punto más distante hidráulicamente; (b) tiempo 
desde el final de la lluvia en exceso hasta que 
finaliza el escurrimiento directo. 

Grimaldi et al. (2012) indican que lo escaso 
de las mediciones de lluvia y su respuesta en 
escurrimiento directo ha vuelto poco práctica la 
estimación del Tc, incluso en cuencas pequeñas, 
y el uso de las ocho definiciones del Tc ha con- 
ducido a estimaciones no consistentes con un 
único método de cálculo. Por ello, en la práctica, 
el Tc se estima haciendo uso de fórmulas empíri- 
cas, que si no son aplicadas teniendo en cuenta 
de manera estricta su intervalo de aplicación 
conducirán sin duda a evaluaciones erróneas. 


Ecuaciones empíricas del Te más utilizadas en 
cuencas rurales 


Temez (1978) expone uno de los primeros con- 
trastes, al menos en lengua hispana, de varias 
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fórmulas del Tc, T, y Tp usando datos reales de 
29 cuencas usadas por el US Army Corps of 
Engineers para deducir su ecuación del 7. En- 
cuentra que las ecuaciones contrastadas tienen 
la expresión general siguiente: 


Tua-| E ] a) 


ser 


siendo L la longitud del cauce principal en 
km; Sc, su pendiente promedio, adimensional. 
Temez (1978), con base en tales datos reales, 
obtiene una ecuación del Tc, cuyos parámetros 
son a: = 0.30, ff = 0.50 y y = 0.76. 

Grimaldi ef al. (2012) contrastan las seis fór- 
mulas siguientes: Giandotti, Johnstone-Cross, 
Kirpich, California Highways and Public Works, 
Natural Resources Conservation Service y Vi- 
parelli, esta última igual al cociente de L entre 
una velocidad media del flujo, que varía de 3.6 
a 5.4 km/h. El contraste se realiza en cuatro 
cuencas del estado de Texas, EUA, con áreas 
de cuenca variando de 13 a 120 km, y 47, 19, 
52 y 20 eventos de lluvia-escurrimiento en cada 
una; encuentran una gran dispersión en los Tc 
calculados en cada evento. Temez (1978) y Gri- 
maldi et al. (2012) encuentran que la fórmula de 
Giandotti sobreestima el Tc en cuencas grandes. 

Respecto a la ecuación de Kirpich, Grimaldi 
et al. (2012) encuentran, al igual que Gericke 
y Smithers (2014), que subestima el Tc en las 
cuencas chicas. De acuerdo con Fang, Thomp- 
son, Cleveland, Pradhan y Malla (2008), Kirpich 
desarrolló su ecuación con datos de cuencas 
pequeñas del estado de Tennesse, EUA, que va- 
riaron de 4 hectáreas a 0.45 km?, y tuvieron pen- 
dientes del 3 al 12%; señalan que sus resultados 
tiene sesgo en cuencas con flujo preponderante 
en cauces. Williams (1950) indica que el cociente 

L INSc , con Sc estimada como el desnivel total 
del cauce (H) en metros entre L en km, no debe 
exceder de 305 para estar dentro de intervalo de 
los datos reales usados por Kirpich. 

En resumen, el polemista sugiere que si 
el autor desea conservar al Te como variable 
predictiva, que se integre una relación de sus 
fórmulas empíricas y se apliquen con estricto 
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apego al intervalo que tuvieron los datos reales 
utilizados en su deducción para obtener diver- 
sas estimaciones del Tc, y entonces adoptar un 
valor probable y /o confiable. Este planteamien- 
to ha sido indicado y expuesto por Campos- 
Aranda (2015). Gericke y Smithers (2014) 
presentan estas relaciones o listas exhaustivas 
para el Tc, T, y Tp. 
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RESPUESTA AL POLEMISTA 


El artículo se centra en dos aspectos: 1) cómo 
podemos estimar los gastos máximos anuales 
asociados con distintos periodos de retorno en una 
cuenca para la que se cuenta con mediciones 
en su estación hidrométrica, y 2) cómo pode- 
mos hacerlo en sitios en donde no se cuenta 
con mediciones. Para ambos casos se utiliza 
como parámetro base el promedio de los gastos 
máximos anuales registrados en las estaciones 
hidrométricas. 

En el primer caso, al dividir cada valor de los 
gastos máximos anuales registrados en una es- 
tación hidrométrica entre su media se obtienen 
valores transformados, cuya media es 1 en todas 
las estaciones; de tal manera que si se conside- 
ran aquellas en las que la variabilidad relativa 
a dicha media (el coeficiente de variación) es 
similar, se obtendrán muestras con la misma 
media y con un segundo momento semejante. 
Al pensar en los aspectos que inciden en esa 
variabilidad relativa, consideramos que, por un 
lado, en general, conviene no agrupar cuencas 
(estaciones) de diferente región hidrológica 
(RH) y, por otro, dentro de cada RH, analizar 
los aspectos que pueden incidir en la variabi- 
lidad (en relación con su media) de los gastos 
máximos registrados en las distintas estaciones 
(p. ej., la forma en que inciden los fenómenos 
meteorológicos y muy particularmente si son 
afectados de manera directa o no por los hu- 
racanes). 

Dado que la hipótesis implícita es que las 
muestras ampliadas provienen de la misma 
población (se consideraron como opciones las 
funciones Gumbel y Doble Gumbel), se verificó 
la hipótesis primero con la prueba de Fisher; 
cuando dicha prueba llevaba a rechazarla, se 
generaron 10 registros sintéticos del tamaño 
de la muestra histórica, verificando que se ob- 
tuviera un rango de coeficientes de variación, 
que incluyera el rango correspondiente a las 
muestras históricas. Queremos enfatizar que no 
hablamos de “homogeneidad hidrológica”, que 
consideramos un concepto demasiado amplio e 
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impreciso, y que sólo nos interesa la “homoge- 
neidad” en relación con las variaciones de los 
máximos de cada año respecto a su media. 

En cuanto al punto 2), utilizamos nuevamen- 
te como parámetro a estimar a la media de los 
gastos máximos anuales, y consideramos que 
dicha media depende de indicadores de lo que 
llueve en la cuenca (su área y su precipitación 
o el producto de ambos), de las pérdidas y de 
la velocidad de respuesta de la cuenca. Los 
indicadores adecuados de las pérdidas y de la 
velocidad de respuesta no son evidentes; en 
particular, los indicadores de la velocidad de 
respuesta (la longitud del cauce principal, el 
tiempo de concentración, por ejemplo) están 
correlacionados con el área de la cuenca, de 
tal manera que, efectivamente, en la ecuación 
global de ajuste existe una correlación espuria. 
Se recomienda por ello que al seleccionar la fun- 
ción de ajuste se verifique que los exponentes de 
cada indicador tengan un sentido lógico. 

Queremos resaltar que el trabajo que se 
está realizando ha implicado un arduo proceso 
de depuración de la información, que no sólo 
incluye la verificación de la independencia y 
homogeneidad de las series sino varias cosas 
más, dentro de las que destaca la verificación de 
valores muy grandes o muy chicos (en el primer 
caso, si ocurrió una tormenta importante, si en 
las cuencas vecinas también se registraron va- 
lores importantes, etc., y en el segundo de una 
forma más cualitativa). También se verificó que 
las cuencas asociadas con las estaciones hidro- 
métricas no tuvieran Obras importantes aguas 
arriba, es decir, que se conservara la naturaleza 
aleatoria de los escurrimientos. 


Comentarios específicos (procurando 
seguir el orden en que expone el polemista) 


Efectivamente, no se expone en el artículo cómo 
se verificó la homogeneidad y la independencia 
de los registros, lo que no significa que no se 
hiciera un análisis al respecto. En el primer caso, 
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se detectaron, por ejemplo, los casos de las cuen- 
cas del oriente del valle de México en los que las 
actividades del Plan Lago de Texcoco han hecho 
que los gastos máximos anuales se reduzcan de 
forma considerable; respecto a la independen- 
cia, es claro que al tratarse de registros de gastos 
máximos anuales, éstos son independientes 
(excepto si, por ejemplo, el máximo de un año 
ocurre a fines de diciembre y el del año siguiente 
a principios de enero). Cabe comentar que, en 
nuestra opinión, son ya muchos los trabajos que 
ocupan un gran espacio describiendo cómo a 
cada registro se le aplican varias pruebas de 
homogeneidad y de independencia. 

Estamos de acuerdo en que el inciso de la 
página 18 titulado “Agrupación de regiones”, 
se debió llamar “Agrupación de estaciones 
hidrométricas”. En cambio, no estamos de 
acuerdo en que no se deben formar grupos 
dentro de las regiones; de hecho, si no se hace 
esa subdivisión, la hipótesis de que los valores 
transformados provienen de la misma población 
no se podría sostener. 

Respecto a la sugerencia de tomar en cuenta 
una verificación de la homogeneidad hidroló- 
gica regional con base en la fecha de ocurrencia 
media y la regularidad estacional de los gastos 
máximos anuales, consideraremos estudiarla, 
aunque a priori no vemos su sentido. 

Como ya se mencionó, se utilizaron como 
opciones sólo las funciones Gumbel y Doble 
Gumbel; la elección de una de ellas se realizó 
por inspección visual de la gráfica que relaciona 


x = —Ln(In(Tr/(Tr-1) con los gastos máximos. 
Cuando en esa gráfica se observa la conforma- 
ción de dos poblaciones usamos doble Gumbel 
(además, se verificó que los puntos de la segun- 
da población provinieran de eventos meteoroló- 
gicos extremos, en particular huracanes o lluvias 
extremas de invierno). 

En efecto, la obtención de exponentes po- 
sitivos en las dos variables predictivas S y Te 
implica una relación incorrecta de la ecuación 
empírica de la media del Qma. Así lo decimos en 
las conclusiones y recomendaciones. 

El polemista tiene razón en que faltó una 
explicación respecto a las figuras 3 y 5; en la 
primera se habían considerado 15 estaciones 
y se consignó el coeficiente de determinación 
que da el EXCELO directamente al ajustar 
una función potencial (que es mayor al que se 
obtiene en forma más estricta con la ecuación 
R? = (variancia de y - variancia del error)/ va- 
riancia de y. 

La figura 5 tiene menos puntos debido a que 
se eliminaron algunas estaciones hidrométricas 
de las cuencas de la cuenca del Cutzamala por 
el manejo artificial de esa cuenca, así como la 
estación hidrométrica Mezcala, cuyos datos son 
prácticamente los mismos que los de la estación 
Santo Tomás. 

Finalmente, agradecemos al polemista la 
disertación sobre el tiempo de concentración y 
cómo estimarlo; para nosotros sólo es un indica- 
dor que involucra en un solo parámetro (1/s%?) la 
velocidad de respuesta de la cuenca. 
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